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by 

Jack  K.  Hale 


Abstract 


Recent  results  on  the  stability  of  equilibrium  solutions 
of  a  parabolic  equation  are  given  with  indications  of  the 
proofs.  Particular  attention  is  devoted  to  dependence  of 
the  stability  properties  on  the  shape  of  the  domain  and 
the  manner  in  which  nonhomo geneous  stable  equilibria  can 
occur  through  a  bifurcation  induced  by  varying  the  domain. 
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In  this  paper,  we  present  a  few  recent  results  on  the 
asymptotic  behavior  of  the  solutions  of  a  parabolic  equation  of 
the  form 


u^  =  Au  +  f(u)  in  n 

=  0  on  an 


(1) 


where  is  a  bounded  open  set  in  IR^  with  smooth  boundary. 

We  also  discuss  how  the  qualitative  behavior  of  the  stable 
equilibria  depend  upon  the  shape  of  and  the  nonlinear  function 
f.  The  function  f  is  supposed  to  satisfy  conditions  which 
ensure  that  Eq.  (1)  defines  a  strongly  continuous  semigroup 
T^(t)  on 

We  remark  that  the  boundary  conditions  in  (1)  are  not 
important  as  far  as  the  spirit  of  the  problems  to  be  discussed. 

Of  course,  the  details  will  depend  generally  in  a  very  significant 
manner  upon  the  boundary  conditions. 

Let  E£(0)  be  the  set  of  equilibrium  solutions  of  (1); 
that  is,  the  set  of  solutions  of  the  equation 


Au  +  f(u)  *0  in  n 


3u  _ 
5n  " 


0 


(2) 


on  dn  . 


2. 


If  Jl  =  (0,X)  c]R,  then  the  set  E^CO^X)  is  given  by 

EfCO.X)  =  {periodic  solutions  of  period  2X  of  the  (3) 
equation  u^^  +  f(u)  =0} 

In  the  following,  we  let  W®((J))  denote,  res- 

spectively,  the  unstable  and  stable  manifolds  for  an  equilibrium 
point  (|>  of  (1).  The  following  result  for  n  =  1  is  due  to  a 
number  of  authors.  The  references  are  in  the  proof. 

Theorem  1.  (n=l)  If  =  (0,X)  c  ]R,  then 

(i)  the  u-limit  set  of  any  bounded  solution  of  (1)  is 
in  Ej(0,X). 

(ii)  the  u-limit  set  of  any  bounded  solution  of  (1)  is 
a  single  point  in  E£(0,X). 

Ciii)  the  only  stable  equilibrium  points  of  (1)  are  constants. 

fU 

(iv)  If  fCs)ds  •*■•«>  as  u  ±«>,  then  every  solution 

h 

of  (1)  is  bounded.  If,  in  addition,  E£(0,X)  is  a 
bounded  set,  then  there  is  a  maximal  compact  in¬ 
variant  set  A£(0,X)  of  (1), 

Af(0,X) 

A£(0,X)  is  uniformly  asymptotically  stable  and,  for 
any  bounded  set  B  in  H^(n),  dist(T£(t)B,A£(0,X))  -►  0 
as  t  -*•  eo. 


3. 


(v)  If,  in  addition  to  the  hypothesis  in  (iv),  all 

(|»  €  E^COtX)  are  hyperbolic,  then,  for  any  bounded 
set  B  in  the  set 

is  open  and  dense  in  B. 

Proof:  (i)  This  is  due  to  Chafee  [4]  and  is  independent  of  the 

boundary  conditions.  The  idea  is  very  simple. 


If 


V(u) 


f (s)ds)dx 


then  the  derivative  of  V  along  the  solutions  of  (1)  satisfies 

V  (u)  =  -  f  u^  dx  <  0 . 

•'0  ^ 

Since  every  bounded  orbit  is  precompact,  a  simple  application 
of  the  invariance  principle  implies  the  result. 

(ii)  This  result  is  due  to  Matano  [13]  and  is  independent 
of  the  boundary  conditions.  He  used  a  rather  sophisticated 
application  of  the  invariance  principle.  We  sketch  a  proof  based 
on  the  theory  of  dynamical  systems.  The  details  will  appear  in 
Hale  and  Massatt  [6].  The  idea  is  very  simple  and  can  be  traced 
to  Malkin  [12],  Hale  and  Stokes  [7]  and  perhaps  even  further. 


If  (p  is  an  element  of  an  u-lirait  of  an  orbit  which  is  not  a 
single  point,  then  <p  belongs  to  a  continuum  in  E£(0,x)  and 
(|)  is  not  hyperbolic.  The  linear  variational  equation  about  (J) 
has  the  simple  eigenvalue  zero  with  all  other  eigenvalues  being 
nonzero.  Thus,  (|)  belongs  to  a  smooth  one  dimensional  sub¬ 
manifold  M  of  E£.  In  modern  terminology,  (j)  is  normally 
hyperbolic.  One  now  can  show  that  any  solution  of  (1)  which 
remains  in  a  sufficiently  small  neighborhood  of  M  for  t 
sufficiently  large  must  be  of  W^(M),  the  stable  manifold  of 
M.  Finally,  one  shows  that  each  orbit  in  W® (M)  approaches  a 
single  point. 

(iii)  This  result  is  due  to  Chafee  [4] .  The  following 
proof  is  taken  from  a  preprint  of  Bardos,  Matano  and  Smoller  [1]. 
Suppose  (p  is  a  nonconstant  equilibrium  solution  of  (1).  Then 

V  =  d(j)/dx  ;^0,  v  =  0  at  x  =  0,x=X,  and  v  satisfies  the 
equation 

’'xx  *•  ■  '>■ 

» 

Let  Oj^,  Ojj  be  the  spectrum  of  this  differential  operator  with, 
respectively,  homogeneous  Neumann  and  Dirichlet  boundary  conditions 
Since  inf  Oj^  <  inf  Op  and  0  €  Op,  the  result  is  proved. 

(iv)  Using  the  function  V(u)  in  part  (i) ,  one  easily 
shows  that  every  solution  is  bounded.  Thus,  the  w- limit  set  of 
every  solution  of  (1)  belongs  to  E£(0,X).  Since  E£(0,X)  is 


5. 


bounded,  there  is  a  bounded  set  B  such  that  every  solution  o£ 

(1)  eventually  enters  B;  that  is,  Eq,  (1)  is  point  dissipative. 
Since  the  semigroup  T£(t)  is  compact  for  t  >  0,  the  results 
follow  from  the  general  theory  of  dissipative  processes  (see,  for 
example.  Hale  [9,10]). 

(v)  This  result  is  due  to  Henry  [11],  The  idea  of  the 
proof  is  to  observe  first  that  the  solution  operator  for  the 
linear  variational  equation  of  (1)  about  any  point  is  one-to-one. 
This  can  then  be  used  to  show  that,  for  any  (j)  e  E£(0,X)  which 
is  not  stable,  the  set  W^Ccf)  n  B  is  nowhere  dense  in  B  for 
any  bounded  set  B.  As  remarked  by  Mane,  the  one-to-oneness  of 
this  solution  operator  implies  that  the  stable  manifold  W^(({») 
can  actually  be  given  globally  a  manifold  structure.  This  also 
gives  a  proof  of  the  assertion  in  (v) . 

For  n  =  (0,X),  we  have  remarked  that  the  set  E£(0,X) 
coincides  with  the  set  of  2X-periodic  solutions  of  u  +  fCu)  =  ' 
For  any  a  6]R,  let  u(x,a)  be  the  solution  of  this  equation  with 
u(0,a)  =  a,u^(0,a)  =0.  If  a  is  such  that  u(x,a)  is  periodic 
in  X,  let  ZX£(a)  be  the  period.  The  function  u(*,a)  e  E£(0,X 
For  f(u)  an  arbitrary  cubic  polynomial  in  u,  Smoller  and 
Wasserman  [15]  have  shown  that  the  function  X£(a)  has  a  finite 
number  (<2)  of  maxima  and  minima  and  the  second  derivative  of  X 
at  these  points  is  different  from  zero;  that  is,  X£Ca)  is  a 


(a)). 


Morse  function. 


The  above  result  has  important  implications  for  the 
applications.  In  fact,  for  f  any  cubic  polynomial  in  u,  and 
for  X  fixed  and  different  from  the  maxima  and  minima  of  the 
function  X£(a),  the  set  E£(0,x)  consists  only  of  hyperbolic 
points.  For  x  equal  to  one  of  the  extreme  values  of  X£(a), 
there  is  a  bifurcation  of  the  saddle-node  type. 

The  following  qualitative  result  of  Brunovsky  and  Chow  [2] 
has  recently  been  proved. 

Theorem  2.  There  is  a  residual  set  3^  £  <if^QR)  with  the  Whitney 
topology  such  that,  for  any  f  €  3^ ,  the  function  X£(a)  above 
is  a  Morse  function. 

The  proof  is  not  trivial  because  the  function  f  depends 
only  on  u  and  not  on  (x,u).  The  proof  is  based  on  a  detailed 
analysis  of  an  analytic  expression  of  X£(a)  as  a  function  of 
f,a.  It  is  not  a  trivial  exercise  because  there  is  no  simple  way 
to  determine  the  qualitative  properties  of  the  derivatives  of  this 
function  in  a  from  the  derivatives  of  f.  In  fact,  there  are 
nonlinear  functions  f  for  which  X£(a)  is  constant  (see,  for 
example,  Urabe  116]). 

Theorem  2  can  be  appropriately  generalized  to  other  boundary 
conditions  (see  Brunovsky  and  Chow  12]).  Smoller  and  Wasserman 
[15]  have  also  considered  other  boundary  conditions. 

Our  next  objective  is  to  discuss  the  extent  to  which  the 
above  results  are  valid  for  a  bounded  set  0  in  1R”, 


Theorem  5.  If  Ji  is  a  bounded  set  in  ]R“  with  smooth  boundary, 
then 

(i)  the  oj-limit  set  of  a  bounded  orbit  is  in  E£(J2). 

(ii)  the  w-limit  set  of  a  bounded  orbit  is  a  single  point 
if  the  following  condition  is  satisfied: 

If  <ti  €  E£(i^)  is  not  hyperbolic  and  k  is  the 
dimension  of  the  null  space  of  the  operator  A  +  f '  ((j))  in  fi 
with  homogeneous  Neumann  conditions,  then  <}>  belongs  to  a  smooth 
submanifold  of  dimension  k. 

(iii)  convex  implies  the  only  stable  equilibrium  points 
are  constants. 

Civ)  Same  statement  as  (iv)  in  Theorem  1  holds. 

(v)  Same  statement  as  (v)  in  Theorem  1  holds. 

Proof:  The  proof  of  (i) ,  (ii) ,  (iv)  and  (v)  are  essentially  the 
same  as  the  corresponding  assertions  in  Theorem  1. 

(iii)  This  result  was  independently  discovered  by  Casten 
and  Holland  [3],  Matano  [14].  The  proof  exploits 
special  properties  of  the  Laplacian  on  convex  regions 
to  prove  that  the  linear  variational  equation  has  a 
negative  eigenvalue  for  any  nonconstant  equilibrium. 

An  analogue  of  Theorem  2,  as  far  as  hyperbolicity  of 
equilibrium  and  saddle  node  bifurcations,  is  not  known  for  n  in 
IR”  and  seems  to  be  rather  difficult. 


8. 


We  remark  that  part  (iii)  of  Theorem  3  is  also  valid  for 
n  convex  and  the  equations 

=  Au  +  f(u),  in 
=  -g(u,v) 

=  0  on  352 

provided  that  the  spectrum  of  the  operator  9g(C,n)/8v  belongs 
to  the  set  {z:Rez  >0}  uniformly  in  C,n  (see  Bardos,  Matano 
and  Smoller  [1]). 

The  remainder  of  the  discussion  centers  around  the  case 
when  52  is  not  convex  and  the  objective  is  to  understand  more 
about  the  set  of  stable  equilibrium.  Before  doing  this,  we  make 
the  important  remark  that,  when  52  is  convex,  the  qualitative 
structure  of  the  stable  equilibria  is  independent  of  the  non¬ 
linearity  f.  When  52  is  not  convex,  this  will  no  longer  be 
the  case. 

The  following  result  is  due  to  Matano  [M]. 

Theorem  4.  There  is  a  nonlinear  function  f  and  52  cIR”,  n  >  2, 
such  that  (1)  has  a  stable  nonconstant  equilibrium. 

Idea  of  proof:  Suppose  f  has  zeros  only  at  a  <  0  <  b,  they 
ru 

are  simple  and  f(s)ds  -*■-«»  as  u  ±®.  Suppose  52  has  the 

0 
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shape  shown  in  Figure  1  and  let  X2  be  the  minima  of  the 


Figure  1 

second  eigenvalues  of  the  Laplacian  on  and  ^2*  Matano  [  ] 

2 

gives  a  specific  continuous  function  GilR  -»■  IR  such  that  the 
set  where  G  <  0  is  nonempty  with  the  property  that,  if  X2  >^3 
are  such  that  G(X2,  meas  <  0,  then  there  is  a  nonconstant 
stable  equilibrium  of  (1).  For  fixed  be  shows  there  is 

an  such  that  the  above  inequality  is  satisfied. 

The  idea  of  the  proof  is  the  following.  Let  Y  be  the 
subset  of  functions  u  in  such  that  f  u  >  0,  [  u  <  0. 

If  meas  is  small  enough,  it  is  then  shown  that  the  set  Y  has 

a  certain  invariance  property  with  respect  to  T£(t).  A  careful 

application  of  the  maximum  principle  gives  a  minimal  equilibrium 

solution  V  in  Y  stable  from  below  and  a  maximal  equilibrium 

m 

solution  Vj^  in  Y  stable  from  above.  If  it  were  known  that 


there  are  only  a  finite  number  of  equilibrium  solutions,  then  the 
proof  would  be  complete.  Since  this  is  not  known,  another  argument 
must  be  used.  Matano  first  proves  that  any  solution  unstable 
from  above  must  be  strongly  unstable  from  above  in  the  sense 
that  it  can  be  isolated  from  equilibrium  solutions  from  above 
uniformly  in  Y.  He  then  puts  an  ordering  on  the  solutions  from 
above,  uses  Zorn's  lemma  and  the  above  property  of  solutions  unstable 
from  above . 

In  the  proof  of  Theorem  4,  the  nonlinear  function  f  has 
the  three  simple  zeros  at  a  <  0  <  b.  The  equilibrium  points 
a,b  are  stable  and  zero  is  a  saddle  point.  A  stable,  nonconstant 
equilibrium  solution  was  shown  to  exist.  The  argument  of  Matano 
can  be  used  to  show  there  must  be  another  nonconstant  equilibrium 
solution  with  u  <  0,  u  >  0.  Thus,  there  are  at  least 

four  stable  equilibrium  solutions  and  one  unstable  equilibrium 
solution.  This  is  impossible  dynamically  and  there  must  be  some 
other  equilibrium  solutions  which  are  unstable.  In  fact,  an  index 
argument  implies  there  must  be  at  least  three  unstable  equilibria. 
Using  more  of  the  detailed  information  from  the  paper  of  Matano, 
one  can  show  there  must  be  at  least  five  unstable  equilibrium 
solutions.  Consequently,  there  are  at  least  nine  equilibrium 
solutions  for  this  nonlinear  function  f  and  region  fi. 

The  basic  problem  is  to  understand  in  more  detail  how 
variations  in  the  shape  of  the  domain  U  cause  these  additional 
solutions  to  appear.  We  now  summarize  some  work  of  Hale  and  Vegas 
[8]  which  give  a  possible  explanation. 


Let  us  begin  with  an  intuitive  discussion  of  how  the 

stable  nonconstant  equilibrium  solutions  could  appear  as 

secondary  bifurcations.  Suppose  y  e  [0,<»),  5;^  is  a  bounded 

set  in  IR  with  smooth  boundary  with  the  property  that  JJq  is 

convex  and  the  second  eigenvalue  X2(y)  of  -A  on 

is  a  monotone  decreasing  function  of  y,  approaching  zero 

as  y  Also,  suppose  the  third  eigenvalue  of 

-A  on  satisfies  Xj(y)  ^  6  >  0  for  all  y.  Let 

f(v,u)  =  V  u  -  u  ,  V  >  0,  and  let  A  be  the  maximal  compact 

V  >  y 

invariant  set  for  the  equation 

u^  =  Au  +  f(v,u)  in 

t  fJ 

-  0  in  3!!^. 

Fix  V  sufficiently  small  so  that  v  <  6  and  the  only  equilibrium 
solutions  in  Aq  are  the  constant  functions  0,  ±v.  The  set 
A  «  is  then  the  constant  functions  ±v  and  the  unstable  manifold 

V  ^  U 

of  0,  which  is  one  dimensional.  Let  Uq  be  such  that 
XzCUo)  ®  f ' (0)  =  V  .  At  the  point  y  =  Uq,  the  zero  solution 
bifurcates  creating  two  new  equilibrium  solutions  which  are  unstable. 
They  are  unstable  because  the  unstable  manifold  of  zero  becomes  two 
dimensional  -  the  direction  of  bifurcation  is  independent  of  the 
direction  of  the  original  unstable  manifold  in  A  for  y  <  y«. 

V  >y  u 

The  set  A^  ^  for  M  >  Uq  close  to  pg  is  then  two  dimensional 

with  three  unstable  and  two  stable  equilibriums. 
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Now  suppose  that  has  the  shape  shown  in  Figure  1 

and  that  meas  -»■  0  as  p  -v  <».  Then  we  can  find  a 
such  that  the  inequality  in  the  proof  of  Theorem  4  is  satisfied. 
Thus,  there  is  a  stable  nonconstant  equilibrium  solution.  It 
is  conjectured  that  this  occurs  as  a  secondary  bifurcation  from 
the  unstable  nonconstant  equilibria  discussed  above.  In  Figures 
2  and  3  we  have  depicted,  respectively,  the  set  A  as  a 
function  of  p  and  the  conjectured  bifurcation  diagram. 


Figure  2 


Figure  3 
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To  do  things  more  analytically,  we  have  tried  to  discuss 
the  neighborhood  o£  y  =  «  treating  it  as  a  bifurcation  problem 
from  a  double  eigenvalue  zero.  More  specifically,  suppose  y  =  e" 
and  the  region  is  shown  in  Figure  4,  two  circles  ^2 

a  small  channel  between  them. 


Figure  4 

For  e  *=  0,  it  is  clear  there  are  nine  solutions  consisting  of 
all  combinations  of  0,  ±v  on  and  ^2*  Five  of  these  are 

saddles  and  four  are  stable  nodes.  Let  ^2^^^  second 

eigenvalue  of  -A  on  and  let  w^  be  a  unit  eigen¬ 
vector  corresponding  to  Let  u^  be  the  constant  function 

(meas  that  is,  u^  is  a  unit  eigenvector  for  the  eigen¬ 

value  0  of  the  Laplacian  on 

One  can  now  show  that  it  is  possible  to  apply  the  method 
of  Liapunov- Schmidt  for  the  solutions  of  (1)  near  u  =  0  for 
e,v  near  zero.  More  specifically,  for  real  a,&  sufficiently 
small  and  e,v  sufficiently  small,  there  is  a  function  u*(a,3»e,v) 
continuously  differentiable  in  a,3»v  and  continuous  in  e  such 
that  u*(0, 0,0,0)  -  0,  3u*(0,0,0,0)/3(a,3)  -  0  and 
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Au*  +  £(v,au  +6w  +u*)  -  irfCv.au  +3w  +u*)  =  0 

&  &  w  w 

where  iru  is  the  projection  of  u  onto  the  span  of  the  constant 

function  u  and  the  function  w  ;  that  is, 
e  e 

TTU  =  u  +  w  I  w  u. 
e  e 

If  u*(a,6,e,v)  satisfies  the  above,  then 
u  =  au^  +  +  u*(a,3, e,v) 

is  a  solution  of  (1)  if  and  only  if  (a,B»e,v)  satisfy  the  bi¬ 
furcation  equations 

f(v,au  +6W  +u*(a,3,E,v))dx  =  0 

•’fig 

I  w  f(v,au  •►gw  •••u*(a,0,e,v))dx  =  0 

Jo  8  8  8 

8 

are  satisfied. 

If  we  let  (a,g)  =  y»  then  one  can  show  that  these  equations 
have  the  form 

cCy)  *  eLj^Y  ■*^''L2Y  h.o.t.  =  0  (4) 

where  h.o.t.  denotes  higher  order  terms  in  y»e,v, 
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constant  two  by  two  matrices  and  c(y)  is  a  homogeneous  cubic 
two  vector  in  the  components  of  y.  One  can  now  apply  the  method 
in,  for  example,  Chow,  Hale  and  Mallet-Paret  {5]  to  obtain  the 
complete  bifurcation  diagram  for  the  solutions  of  (4),  These 
are  shown  in  Fiture  5.  Figure  5a  shows  the  number  of  solutions 


V 


solution 
) — 


for  a  fixed  (e,v)  and  Figure  5b  shows  the  way  the  solutions 
bifurcation  as  a  function  of  e  for  a  fixed  v. 

In  the  verification  of  the  previous  results,  it  is  crucial 
to  show  that  the  third  eigenvalue  of  the  Laplacian  on  is 

bounded  away  from  zero.  It  would  be  very  interesting  to  obtain 
general  geometric  conditions  on  a  region  to  have  this  latter 

property  satisfied. 
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